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ABSTRACT 

We consider a stationary axisymmetric force- free degenerate magnetosphere of 
a rotating Kerr black hole surrounded by a thin Keplerian infinitely-conducting 
accretion disk. We focus on the closed-field geometry characterize by a direct 
magnetic coupling between the disk and the hole's event horizon. We first argue 
that the hole's rotation necessarily limits the radial extent of the force-free link 
on the disk surface: the faster the hole rotates, the smaller the magnetically- 
connected inner part of the disk has to be. We then show that this is indeed the 
case by solving numerically the Grad-Shafranov equation — the main differential 
equation describing the structure of the magnetosphere. An important element 
in our approach is the use of the regularity condition at the inner light cylinder 
to fix the poloidal current as a function of the poloidal magnetic flux. As an 
outcome of our computations, we are able to chart out the maximum allowable 
size of the portion of the disk that is magnetically connected to the hole as a 
function of the black hole spin. We also calculate the angular momentum and 
energy transfer between the hole and the disk that takes place via the direct 
magnetic link. We find that both of these quantities grow rapidly and that their 
deposition becomes highly concentrated near the inner edge of the disk as the 
black hole spin is increased. 

Subject headings: black hole physics — MHD — accretion, accretion disks - 
magnetic fields — galaxies: active 
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1. Introduction 

This paper is devoted to the subject of magnetic interaction between a rotating black hole and 
an accretion disk around it — a topic that has enjoyed a lot of attention among researchers 
in recent years. Magnetic fields are believed to play an important role in the dynamics 
of accreting black hole systems (e.g., Begelman, Blandford, & Rees 1984; Krolik 1999b; 
Punsly 2001). In particular, they can be very effective in transporting angular momentum 
and the associated rotational energy of either the hole or the disk. 

Where and how this transport takes place and to what observational consequences it can 
lead, is partly determined by the global geometry of the magnetic field lines. Conceptually, 
one can think of two basic types of geometry. The first type is the open-field configuration 
shown schematically in Figure 1. The main topological feature here is that there is no direct 
magnetic link between the hole and the disk. Instead, all the field lines are open and extend 
to infinity. Historically, this configuration was the first to have been considered, and it 
has been studied very extensively during the past three decades (see, e.g., Lovelace 1976; 
Blandford 1976; Blandford & Znajek 1977, hereafter BZ77; MacDonald & Thorne 1982, 
hereafter MT82; Phinney 1983; Macdonald 1984; Thorne, Price, & Macdonald 1986; Punsly 
1989, 2001, 2003, 2004; Punsly & Coroniti 1990; Beskin & Par'ev 1993; Beskin 1997; Ghosh 
& Abramowicz 1997; Beskin & Kuznetsova 2000; Komissarov 2001, 2002b, 2004a). The 
reason for this popularity is that this configuration is related to the famous Blandford- 
Znajek mechanism (BZ77) now widely regarded as the primary process powering jets in 
active galactic nuclei (AGN) and micro-quasars. As Blandford and Znajek showed, a large- 
scale, ordered open magnetic field can extract the rotational energy from a spinning black 
hole and transport it to large distances via Poynting flux (a similar process works along the 
field lines connected to the disk). 

The second type of magnetic field geometry is the closed-field configuration, shown in Fig- 
ure 2. Although it has been occasionally discussed in the literature before the last decade 
(e.g., Zeldovich & Schwartzman, quoted in Thorne 1974; MT82; Thorne et al. 1986; Nitta, 
Takahashi, & Tomimatsu 1991; Hirotani 1999), it is only in the last five years that it has at- 
tracted serious scientific attention (e.g., Blandford 1999, 2000, 2002; Gruzinov 1999; van Put- 
ten 1999; van Putten & Levinson 2003; Li 2000, 2001, 2002a, 2002b, 2004; Wang, Xiao, & Lei 
2002; Wang, Lei, & Ma 2003a; Wang et al. 2003b, 2004). The basic topological structure of 
magnetic field in this configuration is very different from that of the open-field configuration. 
The field lines are closed and directly connect the black hole to the disk. In this so-called 
Magnetically-Coupled configuration (Wang et al. 2002), the energy and angular momentum 
are not taken away to infinity, but instead are exchanged between the hole and the disk 
by the magnetic field. Therefore, magnetic coupling, together with the accretion process, 
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controls the spin evolution and the spin equilibrium of the black hole (Wang et al. 2002, 
2003a). In addition, the rotational energy of the hole can be magnetically extracted (just 
like in the Blandford-Znajek process) and deposited onto the disk leading to a change in 
the disk energy-dissipation profile and hence its observable spectral characteristics (Gammie 
1999; Li 2000, 2001, 2002a, 2002b, 2004; Wang et al. 2003a,b). Finally, if the rotating field 
lines are strongly twisted and become unstable to a non-axisymmetric kink-like instability, 
a strong variability of the energy release may result, which would be a possible explanation 
for quasi-periodic oscillations (QPOs) in micro-quasar systems (e.g., Gruzinov 1999; Wang 
et al. 2004). All these phenomena make the closed- field configuration astrophysically very 
interesting. 

Most of the work that has been done on studying the magnetic field structure around an 
accreting Kerr black hole, including the seminal paper by Blandford & Znajek (BZ77), has 
been performed under the assumption that the magnetosphere above the thin disk is ideally 
conducting and force-free. Then, if one also assumes that the system is stationary and 
axisymmetric, the magnetic field is governed by the general-relativistic version of the force- 
free Grad-Shafranov equation (e.g., MT82; for the full-MHD generalization of this equation 
see Nitta et al. 1991; Beskin & Par'ev 1993; Beskin 1997). Since this is a rather nontrivial 
nonlinear partial differential equation (PDE) with singular surfaces and free functions, it is 
generally not tractable analytically, except in some special simple cases, such as the slow- 
rotation limit (BZ77). However, over the past 20 years, a number of force- free solutions for 
the magnetosphere have been obtained numerically, either by solving the Grad-Shafranov 
equation itself (MacDonald 1984; Fendt 1997) or as an asymptotic steady state of force-free 
degenerate electrodynamics (FFDE) evolution (Komissarov 2001, 2002b, 2004a). Until now, 
most of these studies have been done in the context of the open-field configuration, primarily 
because of its relevance to the jet problem. 

In contrast, most of the work on closed-field configurations has been limited to analytic and 
semi-analytic studies of the effects that magnetic link has on the disk radiation profile and 
on the spin evolution of the black hole. The structure of the magnetosphere has not in fact 
been computed self-consistently. These studies have just assumed the existence of the link 
and made some simplified assumptions about the field distribution on the horizon. 

The only exception to this deficiency is the recent work by Uzdensky (2004) where a force-free 
magnetosphere linking a Keplerian disk to a Schwarzschild black hole has been numerically 
computed for the first time. 

In the present paper, we make the next logical step by extending this previous work to the 
more general case of a rapidly rotating Kerr black hole. This is indeed the most important 
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case, not only because real astrophysical black holes are believed to be rotating, but also 
because the nonlinear terms in the Grad-Shafranov equation, especially the toroidal field 
pressure, become large in this case. As a result, even the existence of closed-field solutions 
is not guaranteed. And indeed, one of the main goals of our present study is to determine 
the conditions for existence of such solutions in Kerr geometry. In other words, we aim at 
determining the limitations that the rotation of the black hole imposes on the direct magnetic 
link between the hole and the disk. In addition, by computing the global magnetic field 
structure, we will be able to study the effect of the black hole rotation on the magnetic field 
distribution on the horizon, the poloidal electric current as a function of poloidal magnetic 
flux, and the location of the inner light cylinder, as well as such astrophysically-important 
processes as angular momentum and energy transfer between the hole and the disk. 

In order to achieve the goal of obtaining numerical solutions of the force-free Grad-Shafranov 
equation in Kerr goemetry, we first analize the mathematical structure of this equation. 
In particular, we pay special attention to its singular surfaces (the event horizon and the 
light cylinder) and the corresponding regularity conditions. Thus, we use the light- cylinder 
regularity condition to determine the poloidal current as a function of poloidal magnetic 
flux, similar to way it was done by Contopoulos, Kazanas, & Fendt (1999) for the case of the 
pulsar magnetosphere (see also Beskin & Kuznetsova 2000; Uzdensky 2003; Uzdensky 2004). 
The event-horizon regularity condition, also known as Znajek's (1977) horizon boundary 
condition, is then used to determine the poloidal flux distribution on the horizon. Thus, 
one does not have the freedom to arbitrarily specify any extra boundary conditions at the 
horizon, and hence there is no problem with causality, in line with the reasoning presented 
by Beskin & Kuznetsova (2000) and by Komissarov (2002b, 2004a) (see also Levinson 2004). 

Finally, although in this paper we deal exclusively with large-scale, ordered magnetic fields, 
we acknowledge the difficulty in justifying the existence of such fields around accreting black 
holes (e.g., Livio, Ogilvie, & Pringle 1999), especially in the closed-field configuration. Also, 
as recent numerical simulations (e.g., Hawley & Krolik 2001; Hirose et al. 2004), there may 
be a significant deposition of energy and angular momentum at the inner edge of the disk 
due to small-scale, intermittent magnetic fields connecting the disk to the plunging region 
(see also Krolik 1999a; Agol & Krolik 2000). 

The paper is organized as follows. § 2 describes the mathematical formalism of force-free ax- 
isymmetric stationary magnetospheres in Kerr geometry. In particular, in § 2.1 we introduce 
the Kerr metric tensor in Boyer-Lindquist coordinates and list several general geometric 
relationships for the future use. In § 2.2 we consider steady-state, axisymmetric, degenerate 
electro-magnetic fields and then discuss the force-free condition and the Grad-Shafranov 
equation. In § 2.3 we consider the black hole's event horizon as a singular surface of this 



- 5 - 



equation and discuss the associated regularity condition, which is also known as Znajek's 
horizon boundary condition. In § 3 we present a simple but robust physical argument that 
demonstrates that a force-free magnetic link between a rotating black hole and the disk can- 
not extend to arbitrarily large distances on the disk, we also argue that the maximal radial 
extent of the magnetic link should scale inversely with the black hole's rotation rate in the 
slow-rotation limit. We confirm these propositions in § 4, where we present our numerical 
solutions of the grad-Shafranov equation. Then, in § 5 we discuss the magnetically-mediated 
angular-momentum and energy exchange between the hole and the disk. We then close by 
summarizing our findings in § 6. 



2. Axisymmetric force-free magnetosphere in Kerr geometry — basic 

equations 

2.1. Kerr geometry — mathematical preliminaries 

In this paper we employ Boyer-Lindquist coordinates (t, r, 9, 0) in Kerr geometry. The metric 
of the four-dimensional space-time can be written in these coordinates as 

ds 2 = {g^uj 2 - a 2 )dt 2 - 2ujg H d<t)dt + g„dr 2 + g ee dQ 2 + g^dfy 2 , (1) 

with the components of the metric tensor given by 

« = |VA, (2) 
2aMr 

<* = ~^r, (3) 



s 2 

2 

9rr = ^, 9ee = P 2 , 94,4, = -co 2 , (4) 



where 



p 2 = r 2 + a 2 cos 2 9, (5) 

A = r 2 + a 2 -2Mr, (6) 

S 2 = (r 2 + o 2 ) 2 -o 2 Asin 2 ^, (7) 
E 

w ee — sin 6 . (8) 
P 



Here, M and a G [0; M] are the mass and the spin parameter (specific angular momentum) 
of the central black hole, respectively. (Throughout this paper we use geometric units, i.e., 
we set both the gravitational constant G and the speed of light c to 1). 
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In order to describe the electromagnetic processes around a black hole, we use the 3+1 split 
of the laws of electrodynamics introduced by MT82 (see also Thorne et al. 1986). In this 
formalism, one splits the four-dimensional spacetime into the global time t and the absolute 
three-dimensional curved space, the geometry of which is described by a three-dimensional 
(3D) metric tensor with components given by equation (4). The electromagnetic field is 
represented by the electric and magnetic field 3-vectors E and B measured by local zero- 
angular-momentum observers (ZAMOs; see Thorne et al. 1986). In order to describe these 
vectors, we will use both the coordinate basis {<%} = {e^} and the orthonormal basis {e^} 
[where the Roman index i runs through the three spatial coordinates (r, 9, (/>)]. Because the 
spatial 3D metric tensor g^ is diagonal, these two bases are related via 

ei = v / <fee;, i = r,9,(f) (9) 

(note: there is no summation over % in this expression!). In particular, in the Boyer-Lindquist 
coordinates in Kerr geometry, we have 

" - ' ~r, e = pe<?, e = G7eT. (10) 



We shall also need the following mathematical expressions: the 3-gradient of a scalar function 
/(x) = f(r,9,4>) in the Boyer-Lindquist coordinates is 



- (d r f) e, + - (d e f) e § + - (d+f) , (11) 
p p w r 

and its square is 

IWI 2 = £^W) 2 

i 

= - 2 {drf? + - 2 mr + — 2 m? ■ w 

p z p z w z 

Finally, the 3-divergence of a 3-vector A can be written as 

V • A = = A* + A\\n Vbf) = -1= (VW\^) . , (13) 



\9\ 



where g is the determinant of the 3-D metric tensor: 



p 2 w pS sin 9 , . 

y\ = -i^ = -^^- ( 14 ) 
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2.2. Stationary axisymmetric ideal force-free magnetosphere in Kerr geometry 

As mentioned above, in the 3+1 split formalism of MT82 a magnetosphere of a rotating Kerr 
black hole is described in terms of two spatial vector fields, E and B. Under the assumptions 
that the magnetosphere is: (1) stationary (d t = 0), (2) axisymmetric (d^ = 0), and (3) 
ideally-conducting, or degenerate (E • B = 0), these two vector fields can be expressed in 
terms of three scalar functions, ^(r, 9), £l F (r, 9), and I(r, 9): 

B(r,9) =B pol + B tor , (15) 

where 



and 



where 



B i = v^xV0 = — ^ e er -^ r e § , (16) 

wp wp 

Btor = B~ % = -L-e~ (17) 

v v aw v 



5Vt 

E(r,0) = E pol = -— W, ^ = 0, (18) 



sn = n F -u. (19) 



Here, \l/(r, 9) is the poloidal magnetic flux function, Q F = fl F (^) is the angular velocity 
of the magnetic field lines, and I(r,9) is (2/c) times the poloidal electric current flowing 
through the circular loop r = const, 9 = const. [Note that our definitions of \l/ and / differ 
from the ones adopted by MT82: \1> = ip MT82 /2ir, I = — (2/c)/ M T82-] 

Next, in this work we are interested in the case of a force-free magnetosphere, i.e., a mag- 
netosphere that is so tenuous that the electromagnetic forces completely dominate over all 
others, including gravitational, pressure, and inertial forces. Even though this framework 
has been widely accepted as a primary tool in describing magnetospheres of black holes and 
radio-pulsars, its usefulness and validity near the event horizon has been seriously challenged 
by Punsly (2001, 2003). However, according to the recent MHD simulations by Komissarov 
(2004b), these worries seem to be unfounded. Therefore, we shall still employ the force- free 
approach in this paper. Correspondingly, we shall write the force-balance equation (in the 
ZAMO reference frame) as 

Pe E+jxB = 0, (20) 

where the ZAMO-measured electric charge density p e and electric current density j are related 
to E and B via Maxwell's equations (see MT82). 
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The toroidal component of the force-free equation immediately leads to 

I(r,9) = /(*), (21) 
i.e., the poloidal electric current does not cross poloidal flux surfaces. 

The poloidal component of equation (20), upon using expressions (15)-(18), yields the so- 
called generally-relativistic force-free Grad-Shafranov equation — the main equation that 
governs the system. In this paper we shall use as a starting point the form of this equation 
given in MT82 (i.e., eq. [6.4] of MT82 slightly modified to account for the change in the 
definition of \l/): 



( a r , oil w .„ _ \ 
— 1 — V* + 

\zu z cr / 

a d^> aw 2 



— ^f( V *) + T=* n W = °- ( 22 ) 



This is a nonlinear 2nd-order elliptic partial differential equation (PDE); it determines \I/(r, 9) 
provided that Qf(^) and are known. We can rewrite this equation as follows: 

LHS = aw 2 V ■ ( 5 (a 2 - 6Q 2 w 2 Wv) = RHS = -JJ'(tf) - Snn'JV) ro 2 (W) 2 , (23) 

where a prime denotes the derivative with respect to e.g., = dl/d^f. 

Upon introducing the quantities 

D = a 2 — 5tt 2 zu 2 , (24) 

and 

QM) S ^ = 4 = -^, (25) 
azu 2 zuA A sin 9 

and upon using identity (13), the left-hand side (LHS) of this equation can be written in a 
compact and convenient form 

LHS = Q- 1 [QD(W) l }, = WL+DTOlnQ. (26) 



Using expression (11), we get the Grad-Shafranov equation in the following final form: 

LHS = d r (^^ r ^d e (^^e)+^(^r^d r Q + ^ e deQ) 

= RHS = -JJ'(tt) - 5QQ' F {^) ro 2 (W) 2 . (27) 
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2.3. Regularity condition at the event horizon 

From the Grad-Shafranov equation in the form (27) it is easy to see that, in general, this 
equation has two types of singular surfaces. One of them is the so-called light cylinder (often 
called in the literature the velocity-of-light surface or simply the light surface) defined as a 
surface where D = 0. We shall discuss it in more detail later (see § 4.2). 

There is also another singular surface of the Grad-Shafranov equation: the event horizon 
defined as the surface where 

A = = a . (28) 
This surface will be the main focus of this section. 

As can be seen from equation (6), the event horizon is a constant-r surface, 

r(0) =r H = M + VM 2 - a 2 = const . (29) 

In addition, the frame-dragging frequency u defined by equation (3) is also constant on the 
horizon, 

uj{r = r H , 6) = Vt H = — — = const . (30) 

2Mr H 

This constant is what is conventionally called the rotation rate of the Kerr black hole. 

Because the horizon is surface of constant r, one can immediately see that it is a singular 
surface of equation (27). This is because the coefficient in front of the 2nd-order derivative in 
the direction normal to this surface (in this case, radial) vanishes, even though the coefficient 
in front of the 2nd derivative in the ^-direction does not. 

The fact that the event horizon is just a singular surface of the Grad-Shafranov equation is 
extremely important. It means that one cannot impose an independent boundary condition 
for the function \I/(r, 0) at the horizon. One can only impose a regularity condition there 
(Beskin 1997; Komissarov 2002b, 2004a; Uzdensky 2004). Mathematically, this condition 
means that there should be no logarithmic terms in the asymptotic expansion of \I/(r, 9) near 
r — Th (see MT82). Physically, the regularity condition originates from the requirement 
that freely-falling observers measure finite electric and magnetic fields near the horizon (see 
Thorne et al. 1986). Alternatively, the event horizon regularity condition can be obtained 
from the fast-magnetosonic critical condition in the limit in which plasma density goes to 
zero and hence the inner fast magnetosonic surface approaches the horizon (Beskin 1997; 
Beskin & Kuznetsova 2000; Komissarov 2004a). In the present paper, we will not repeat the 
rigorous derivation of this condition (we refer the reader to MT82 or Thorne et al. 1986). 
Instead, we just note that as a result of the regularity requirement, one expects both the 1st 
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and 2nd radial derivatives of \1/ to remain finite at the horizon. Therefore, when applying 
the Grad-Shafranov equation (27) at r = r H) one can just simply set A = 0. Then, after 
some algebra, one gets: 

where 

tt o (0) = *(r = r H ,0) . (32) 



I 2 IM8)] = fa- ^f) 2 + const , r = r H , (31) 



In the absence of a finite line-current along the axis = 0, i.e., when 1(6 = 0) = 0, the 
integration constant is zero and hence 

/ = ±5tt — , r = r H . (33) 

p da 

As for the choice of sign in this expression, it can be shown that the correct sign must be 
plus [remember that MT82 have minus sign because we define I(^f) with an opposite sign]; 
this comes from the requirement that Poynting flux measured by a ZAMO in the vicinity of 
the horizon is directed towards the black hole (e.g., Znajek 1977, 1978; BZ77; MT82). Thus, 
we have 

7-hT, rav rn rod *o 2Mvh sin 6 d^o 

nMe)] = sn- — = p2 an — , r = r H . (34) 

Equation (34) was first derived by Znajek (1977) and is frequently referred to as the "Zna- 
jek's horizon boundary condition". We stress, however, that, because the event horizon is 
a singular surface of the Grad-Shafranov equation, one cannot really impose a boundary 
condition there; expression (34) actually follows from the Grad-Shafranov equation itself 
under the condition that the solution be regular at r = ra- 
il is interesting to note that, because not only the 2nd-, but also the lst-order radial deriva- 
tives of iff drop out of the Grad-Shafranov equation when A is set to zero, this equation 
becomes an ordinary (as opposed to a partial) differential equation at the horizon! This 
implies that the horizon poloidal magnetic flux distribution, if> (6), is connected to the mag- 
netosphere outside the horizon only through the functions I (iff) and f2p(\I/) and not through 
any radial derivatives. From the practical point of view, this fact means that equation (34) 
can be viewed as a Dirichlet-type boundary condition that determines the function iffo(6) 
once both I (iff) and flp(^f) are given. It is important to emphasize that we really have only 
one relationship on the horizon — equation (34) — between three functions [$ (9), and 
and hence one needs to find some other conditions, set somewhere else, to fix I (iff) 
and £|p(\l/) if one wants to use (34) to calculate ^ffo(d)- We shall return to this important 
point in § 4.1. 
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3. Disruption of the hole disk magnetic link by the black hole rotation 

The main topic of this paper is a force-free magnetic link between a Kerr black hole and 
a thin, infinitely conducting Keplerian accretion disk around it. Thus, we are primarily 
interested in the closed-field configuration depicted schematically in Fig. 2. In contrast to 
the open-field configuration, in which all the field lines piercing the event horizon extend to 
infinity, in the closed-field configuration, magnetic field lines connect the black hole to the 
disk, forming a nested structure of toroidal flux surfaces. In this section we will examine 
the conditions under which this configuration can exist and, in particular, will discuss the 
limitations that the rotation of the black hole imposes on the radial extent of the force-free 
magnetic link between the disk and the hole. 

First, we would like to point out that a magnetically-linked black hole-disk system is dra- 
matically different from a magnetically-linked star-disk system in at least one important 
aspect. Indeed, let us examine the system's evolution on the shortest relevant, i.e., rotation, 
timescale. In the case where the central object is a highly-conducting star, such as a neutron 
star or a young star, it turns out that no steady state configuration with the topology similar 
to that presented in Figure 2 is possible. This is because both the disk and the star can 
be regarded (on this short timescale) as perfect conductors, so that the footpoints of the 
field lines that link the two are frozen into their surfaces. Hence, the disk footpoint of a 
given field line rotates with its corresponding Keplerian rotation rate, Q K (r), whereas the 
footpoint of the same field line on the star's surface rotates with the stellar angular veloc- 
ity Q*. Therefore, each field line connecting the star to the disk [with the exception of a 
single line connecting to the disk at the corotation radius r co where f2(r co ) — ^*] is subject to 
a continuous twisting. This twisting results in the generation of toroidal magnetic flux out 
of the poloidal flux, which tends to inflate and even open the magnetospheric flux surfaces 
after only a fraction of one differential star-disk rotation period (e.g., van Ballegooijen 1994; 
Uzdensky et al. 2002; Uzdensky 2002a,b). 

On the other hand, in the case of a black hole being the central object the situation is 
different. The key difference is that, unlike stars, black holes do not have a conducting 
surface. On the contrary, they are actually effectively quite resistive, in the language of the 
Membrane Paradigm (see Znajek 1978; Damour 1978; MacDonald & Suen 1985; Thorne et al. 
1986). The rather large effective resistivity makes it in principle possible for the field lines 
frozen into a rotating conducting disk to slip through the event horizon. This fact makes a 
quest for a stationary closed-field configuration in the black-hole case a reasonable scientific 
task, since it is at least conceivable that such configurations may in principle exist. 

In our previous paper (Uzdensky 2004) we studied exactly this question for the case of a 
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Schwarzschild black hole. We found that a stationary force-free configuration of the type 
depicted in Figure 2 indeed exists in this case. At the same time, however, there is of course 
no guarantee that a similar configuration will exist in the Kerr case. This is because the 
nonlinear terms in the Grad-Shafranov equation that correspond to field-line rotation and 
toroidal field pressure are no longer small in the Kerr case, whereas in the Schwarzschild 
case these terms, although formally finite, were only at a few per cent level. In fact, we can 
make an even stronger statement: even for a slowly-rotating Kerr black hole, a force- free 
configuration in which magnetic field connects the polar region of the horizon to arbitrarily 
large distances on the disk (which is precisely the geometry depicted in Fig. 2) does not 
exist! We shall now present the basic physical argument for why this must be the case. 

Let us suppose that a force-free configuration of Figure 2, where all the field lines attached to 
the disk at all radii thread the event horizon, does indeed exist. First, let us consider the polar 
region of the black hole, r = r H , 9^0. Suppose that near the rotation axis the flux \&o(#) 
behaves as a power law: ^ ~ # 7 (the most natural behavior corresponding to a constant 
poloidal field being ty ~ 9 2 ). Then note that in a configuration under consideration, the field 
lines threading this polar region connect to the disk at some very large radius ro(\l/) ^> th- 
Since the field lines rotate with the Keplerian angular velocity of their footpoints on the 
disk, f2p(\l/) ~ r 3//2 (\I/) — > as ^ — > 0, one finds that, for sufficiently small \1/ [and hence 
sufficiently large ro( x I / )], f^^) becomes much smaller than the black hole rotation rate 
£l H = a/2r H . Now let us look at the event horizon regularity condition (34). For the field 
lines under consideration, we find that sm9d^f /d9 ~ 99 1 ~ 1 ~ \I> and 5il = f2p(\l/) — ~ 
-Sl H = const ^ 0, * -> 0. Thus, 

n H V atf , as -> , (35) 

and, correspondingly, 

//'(* -> 0) ~ a 2 ^ . (36) 

Now, let us look at the force-free balance on the same field lines but far away from the 
black hole, at radial distances of the order of r ~ tq ^> th- At these large distances a ~ 1 
and 5Qw c, so that the electric terms in the Grad-Shafranov equation are small and the 
coefficient D is close to 1. Then the LHS of the Grad-Shafranov equation (27) is essentially 
a linear diffusion- like operator and can be estimated as being of the order of \l//r 2 . We 
see that both the LHS and the RHS given by equation (27) scale linearly with \1> but the 
LHS has an additional factor ~ r~ 2 . Thus we conclude that at sufficiently large distances 
this term becomes negligible when compared with the J/'(\l/)-term (36). In other words, 
the toroidal field, produced in the polar region of the horizon by the black hole dragging 
the field lines along, turns out to be too strong to be confined by the poloidal field tension 
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at large distances. In fact, this argument suggests that the maximal radial extent r max of 
the region on the disk connected to the polar region of a Kerr black hole should scale as 
r max ~ t h/0' i n the limit a — > 0. One should note that, in the Schwarzschild limit a — > 0, 
this maximal distance goes to infinity and hence a fully-closed force-free configuration can 
exist at arbitrarily large distances, in agreement with the conclusions of our paper I. [Also 
note that if one tries to perform a similar analysis for the Schwarzschild case, then from the 
horizon regularity condition one finds that J(tf) = n K (V) sm9((M/d9)\ r=rH ~ • r~ 3/2 (#). 
Then, assuming that r (\P) is a power law at large distances, r ~ r (^), the toroidal- field 
pressure term can be estimated as ~ \I> ■ rjj" 3 (^)- We can thus see that at large 

distances this term becomes negligible compared with the LHS (~ ^r~ 2 ), so no limitation 
on the radial extent of the magnetic link can be derived.] 

We also would like to remark that this finding is not really surprising in view of some im- 
portant properties axisymmetric force-free magnetospheres, known from the general theory 
of the (non-relativistic) Grad-Shafranov equation (see, e.g., van Ballegooijen 1994; Uzden- 
sky 2002b). This analogy is so important that we would like to make a digression to describe 
it here. Let us consider a closed simply-connected (i.e., without magnetic islands) axisym- 
metric configuration like the one shown in Figure 2. Then start to increase gradually the 
overall magnitude (which we shall call A) of the nonlinear source term — the so-called 

generating function — starting from zero. As we are doing this, let's keep the functional 
shape of as well as the boundary conditions for \1>, fixed. Then one finds the following 

interesting behavior: there is a certain maximal value A max (whose exact value depends on 
the details of the functional shape of I(^f) and the boundary conditions), such that one 
finds no solutions of the Grad-Shafranov equation for A > A max . For A < A max , one actually 
finds two solutions and these two solutions correspond to two different values of the field-line 
twist angles A$(\l/). In the limit A <C A max the two solutions are remarkably different. One 
of them corresponds to A$ ~ A/A max <C 1; it is very close to the purely potential closed- 
field configuration and can be obtained as a perturbation from the potential solution. The 
other solution corresponds to some finite distribution A$ c (\l/), in general of order 1 radian, 
and is characterized by very strongly inflated poloidal field lines. This configuration in fact 
approaches the open-field geometry in the limit A — > 0. Now, as one increases A, the dif- 
ference between the two solutions decreases and they in fact merge into one single solution 
at A = A max . The corresponding configuration shows some modest inflation of the poloidal 
field and corresponds to the field line twist angles that are finite (i.e., of order 1 radian) but 
less than A$ c (\l/). Most importantly, as we mentioned above, no solutions with the required 
simple topology (i.e., without magnetic islands) exist for A > A max . 

Clearly, this is exactly what happens in the Kerr black hole case. Indeed, in this case the 
regularity condition (34) requires that the generating function be of the order of 
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a 2 ^ for small In a certain sense, the spin parameter a 2 effectively plays the role of the 
parameter A from our example above. If one considers a configuration in which the magnetic 
link extends to a radius r max on the disk and fixes the disk boundary conditions, it turns out 
that there is a critical maximum value o^ax beyond which no solution can be found. From 
the argument presented in the beginning of this section we expect that a max scale inversely 
with r max ; in particular, for an infinitely extended link (r max — > oo), one finds a max — > and 
no solution is found for any a > 0! 

To sum up, even though the field lines can, to a certain degree, slip through the horizon 
because the latter is essentially resistive, in some situations the horizon is not resistive 
enough to ensure the existence of a steady force-free configuration! Indeed, the field lines are 
"dragged" by the rotating black hole to such a degree that, in order for them to slip through 
the horizon steadily, they must have a certain rather large toroidal component. When, 
for fixed disk boundary conditions, the black-hole spin parameter a is increased beyond a 
certain limit a max (r , inax ), this toroidal field becomes so large that the poloidal field tension is 
no longer able to contain its pressure at large distances. 

Finally, the argument put forward in this section proves that it may be not only interesting 
but also in fact necessary to consider hybrid configurations in which at least a portion of the 
field lines are open and magnetic disk-hole coupling plays a more limited role. 



4. Numerical Simulations 

In order to verify the proposition put forward in the preceding section and to study the 
magnetically-coupled disk-hole magnetosphere, we have performed a series of numerical cal- 
culations. We obtained the solutions of the force-free Grad-Shafranov equation correspond- 
ing to various values of two parameters: the black-hole spin parameter a and the radial 
extent R s of the magnetic link. In this section we shall describe the actual computational 
set-up of the problem, including the boundary conditions and the numerical procedure; we 
shall also present the main results of our calculations. 



4.1. Problem formulation and boundary conditions 

We start by describing the basic problem set-up and the boundary conditions. 

The simplest axisymmetric closed-field configuration one could consider is that shown in 
Figure 2. In this configuration, all magnetic field lines connect the disk and the hole. Fur- 
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thermore, the entire event horizon and the entire disk surface participate in this magnetic 
linkage; in particular, the field lines threading the horizon very close to the axis 9 = are 
anchored at some very large radial distances in the disk: 

^ (9 -> 0) = (r = r H ,9 -> 0) = ^ disk (r -> oo) = *(r -> oo, = tt/2) (37) 

However, as follows from the arguments presented in § 3, a steady-state force- free configu- 
ration of this type can only exist in the case of a Schwarzschild black hole; in the case of a 
Kerr black hole, even a slowly-rotating one (a <C M), such a configuration is not possible. 
And indeed, in complete agreement with this point of view, in our simulations, we were not 
able to obtain a convergent solution even for a Kerr black hole with the spin parameter as 
small as a = 0.05. 

Also in §3 we proposed a conjecture that, for a given value of a, the magnetic link between 
the polar region of the black hole and the disk cannot, generically, extend to distances on 
the disk larger than a certain r max (a). The exact value of r max depends on the details of the 
problem, such as the exact flux distribution \P d (r) on the surface of the disk, etc. However, 
we proposed that r max is a monotonically decreasing function of a, and, more specifically, 
in the limit a — > 0, r max is inversely proportional to a. For a finite ratio a/M = 0(1), we 
expect that the magnetic link can only be sustained over a finite range of radii not much 
larger than the radius of the Innermost Stable Circular Orbit r Isco . 

In order to test these propositions, we set up a series of numerical calculations aimed at 
solving the Grad-Shafranov equation for various values of two parameters: the black-hole 
spin parameter a and the radial extent of the magnetic coupling on the disk surface R s . 

Correspondingly, in order to investigate the dependence on the radial extent of magnetic cou- 
pling, we modified the basic geometry of the configuration by allowing for two topologically- 
distinct regions: region of closed field lines connecting the black hole to the inner part of the 
disk r < R s , and the region of open field lines extending from the outer part of the disk all 
the way to infinity. 2 This configuration is shown in Figure 3. We count the poloidal flux on 
the disk from the radial infinity inward, so that ^(r = oo) = 0, and ^d(f — r n) = ^tot- 
The disk flux distribution may still be the same as in the configuration of Figure 2; however, 
now there is a critical field line ^ s = ^ d (R s ) < ^ tot that acts as a separatrix between open 



2 In general, open field lines originating from the disk may carry a magnetocentrifugal wind (Blandford & 
Payne 1982) and the resulting mass-loading may make a full-MHD treatment necessary for these field lines. 
Here, however, we shall ignore this complication and will assume the force-free approach to be valid in this 
part of the magnetosphere as well. 
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field lines < *ff s ) and closed field lines < ^ < *tot) connecting to the black hole. 
Correspondingly, the poloidal flux on the black hole surface varies from \& = \I> S at the pole 
6 = to * = * tot at the equator 9 = 7r/2. 

It is worth noting that a more general configuration would also have some open field lines 
connecting the polar region of the black hole to infinity. In fact, such a configuration would 
be more physically interesting because these open field lines would enable an additional 
extraction of the black hole's rotational energy via the Blandford-Znajek mechanism (BZ77). 
We shall call this a hybrid configuration because the disk-hole magnetic coupling and the 
Blandford-Znajek mechanism operate simultaneously. In the present paper, however, we 
shall assume that no such hole-infinity open field lines. We make this choice not because 
of any physical reasons but simply because of technical convenience: we want to isolate 
the effect tof disk-hole coupling. In addition, as we discuss in more detail in § 6, a proper 
treatment of these open field lines would require a more complicated numerical procedure 
than that needed for the field lines that connect to the conducting disk. 

In addition to boundary conditions, one has to specify the angular velocity f2p(\l/) of the 
magnetic field lines. Since we assume that the disk is a perfect conductor, and since in our 
field configuration all the field lines go through the disk, this angular velocity is equal to 
that of the matter in the disk. 

Now let us consider the open field lines \P < ^ff s . In principle, since they are attached to a 
rotating Keplerian disk, these lines rotate differentially with the angular velocity = 
^ifjro^)]- Correspondingly, just as the closed field lines going into the black hole or the open 
field lines in a pulsar magnetosphere, they have to cross a light cylinder and therefore have 
to carry poloidal current whose value must be consistent with, and indeed determined 

by, the regularity condition at the light cylinder. Because this outer light cylinder is very 
distinct from the inner light cylinder that is crossed by the closed field lines entering the 
event horizon, we in general would expect the function I(^ff < ^ff s ) be very different from 
the function I (iff > \& s ). In particular, we would expect a discontinuous behavior, J°P en = 
lim I (iff < \tO ^ /closed = jj m > ,|/ \ even though the field-line angular velocity 

Qp = f2fsr[ro(\l/)] remains perfectly continuous and smooth at ^ = ^/ s . 

Dealing with such a discontinuity in I(iff) across the separatrix iff = iff s presents certain 
numerical difficulties, especially taking into account that the location of the separatrix 
r s (9) = r(\& = ^ s ,9) is not known a priori. Therefore, in the present study we decided 
to simplify the problem by introducing the following modifications: we require that the 
outer part of the disk, r > R s , be nonrotating: < if/ s ) = 0. 
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Correspondingly, the open field lines do not cross an outer light cylinder, and so < \I/ S ) = 
0. To put it in other words, we just take the open-field outer part of the disk magnetosphere 
to be potential. Next, in order to avoid the numerically-challenging discontinuities in f2p(\l/) 
and at \& = \I/ S , we slightly modify the disk rotation law just inside of R s by taking Qp 
smoothly to zero over a small (compared with the total amount of closed flux) poloidal flux 
range. In particular, we used the following prescription: 

n F (*) = o, ^<^ s , 

Q F (V) = Q x [ ro (*)]. t an 2 (^^£), vjr > * s , (38) 
where = 0.2(* tot - * s ) and (see equation [5.72] of Krolik 1999, p. 117) 



n * ) = «£hn- (39) 



These modifications enabled us to focus on examining how black hole rotation (i.e., the 
spin parameter a) limits the radial extent R s of the force-free magnetic coupling, while 
at the same time avoiding certain numerical difficulties resulting from the discontinuous 
behavior of poloidal current We believe that these modifications do not lead to any 

significant qualitative change in our conclusions, especially in the case of small a and large 
R s . Nevertheless, we intend in the future to enhance our numerical procedure so that it 
become fully capable of treating this discontinuity. 

Let us now describe the computaional domain and the boundary conditions. 

First, because of the assumed axial symmetry and the symmetry with respect to the equa- 
torial plane, we performed our computations only in one quadrant, described by 6 G [0,7r/2] 
and r G [r# , oo]. Thus, we have four natural boundaries of the domain: the axis 9 = 0, the 
infinity r = oo, the equator 6 = ir/2, and the horizon r = r#. Of these, the axis and the 
equator require boundary conditions for \l>, whereas the horizon and the infinity are actually 
regular singular surfaces and so we only impose regularity conditions on them. 

The boundary condition on the rotation axis is particularly simple: 

yj/(r,e = 0) = ^ s = const. (40) 



The equatorial boundary, 6 = ir/2, actually consists of two parts: the disk (considered to 
be infinitesimally thin) and the plunging region between the disk and the black hole. The 
border between them, i.e., the inner edge of the disk, is assumed to be very sharp and to lie 
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at the ISCO: r in = risco(o); ^isco varies between risco = 6M for a Schwarzschild black hole 
(a = 0) and r ISCO = M for a maximally rotating Kerr black hole (a — > 1). 

Let us first discuss the boundary conditions at the disk surface, r > r in . Depending on the 
resistive properties of the disk, and on the timescale under consideration, one can choose 
between two possibilities, both of which appear to be physically sensible: 

1) If one is interested in time-scales much longer than the characteristic rotation timescale 
but much shorter than both the accretion and the magnetic diffusion timescales, then it is 
reasonable to regard the poloidal flux distribution across the disk to be a fixed prescribed 
function, which must be specified explicitly. Thus, in this case one adopts a Dirichlet-type 
disk boundary condition: 



The function ^(r) is arbitrary; the only requirement that must be imposed in accordance 
with the discussion above is the convention that 4^(r = oo) = and \E'd(r in ) = \& tot . Since 
we don't have any good physical reasons to favor one choice of ^d(r) over any other, we in 
this paper just choose it arbitrarily to be a power-law with the exponent equal to —1: 



2) If one looks for a configuration that is stationary on timescales much longer than the effec- 
tive magnetic diffusion time (while perhaps still much shorter than the accretion time scale), 
then one should regard the disk as effectively very resistive for the purposes of specifying the 
disk boundary condition. This situation may arise in the case of a turbulent disk; for such a 
disk, the effective magnetic diffusivity rj can probably be estimated as i] turh = ass c sh, in the 
spirit of the a-prescription for the effective viscosity in the SS73 model. Then, the character- 
istic radial velocity of the magnetic footpoints across the disk is roughly Vf p ~ assc s (B r / B z )d. 
For the ratio (B r /B z )d of order 1, this velocity is much greater (by a factor of r/h) than 
the characteristic accretion velocity. Therefore, the only way one can have a steady-state 
configuration on the diffusion time-scale (which, according to the above estimate is of the 
order of the disk sound crossing time r/c s ) is for the poloidal field to be nearly perpendicular 
to the disk, B r <C B z . This requirement translates into a simple von-Neumann boundary 
condition for \&(r, 9) at the disk surface: 



*(r > ri„, 9 = tt/2) = * d (r) 



(41) 




(42) 




(43) 



In our present paper, however, we chose the Dirichlet-type boundary condition represented 
by equations (41)-(42) and set ^ t ot — 1 throughout the paper. 



-19- 



In the plunging region (r# < r < r in , 6 = re/ 2) we have chose 



\&(r H < r < r in ,9 = ir/2) = ^ 



tot = 



^d(r in ) = const . 



(44) 



This choice appears to be physically appropriate for an accreting (and not just rotating) 
disk. The reason for this is that the matter in this region falls rapidly onto the black hole 
and thereby stretches the magnetic loops in the radial as well as the azimuthal directions, 
greatly reducing the strength of the vertical field component. The horizontal magnetic field 
then reverses across the plunging region, which is thus described as an infinitesimally thin 
non-force-free current sheet lying along the equator. In essense, this situation is directly 
analogous to the case of a force-free pulsar magnetosphere, where all the field lines crossing 
the outer light cylinder have to be open and extend out to infinity, thus forming an equtorial 
current sheet (Beskin 2003; van Putten & Levinson 2003). In the black-hole case, one could 
still consider an alternative picture of the plunging region with some field lines crossing the 
equator inside r m . However, in this case one would still have to have a non-FFDE equatorial 
current sheet inside the inner light cylinder, as was shown by Komissarov (2002b, 2004a). 

Finally, as we have discussed in § 2.3, the event horizon is a regular singular surface of the 
Grad-Shafranov equation. Correspondingly, one cannot and need not impose an additional 
arbitrary boundary condition here (e.g., Beskin & Kuznetsova 2000; Komissarov 2002b, 
2004a). Instead, one imposes the regularity condition (34); this condition has the form of 
an ordinary differential equation (ODE) that determines the function ^o{@) provided that 
both fip(\l/) and are given. Thus, from the procedural point of view, this condition can 
be used as a Dirichlet boundary condition on the horizon. It is important to acknowledge, 
however, that one does not have the freedom of specifying an arbitrary function ^f(9) and 
then studying how the information contained in this function propagates outward and affects 
the solution away from the horizon. The function ^o(9) is uniquely determined once Of(^) 
and are given and thus there is no causality violation here. 

Similarly, the spatial infinity r = oo is also a regular singular surface of the Grad-Shafranov 
equation and thus can also be described by a regularity condition. In this sense, the horizon 
and the infinity are equivalent (e.g., Punsly & Coroniti 1990). Note that, in our particular 
problem set-up, the sitiuation at infinity is greatly simplified because we have set f2p(\l/) = 
on the open field lines extending from the disk. Because of this, there is no outer light 
cylinder for these lines to cross and thus one can also set J(\l/ < = 0. Then, at very 
large distances (r ^> r#), the Grad-Shafranov equation (27) becomes a very simple linear 
equation: ^/ rr + r~ 2 sin Odei^e/ sin#) = 0, and the asymptotic solution that corresponds to 
the open field geometry with finite magnetic flux is just 



^(r = oo,0) = # s cos6>. 



(45) 
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4.2. Light-Cylinder Regularity Condition 

At this point the problem is almost completely determined. The only thing we still have 
to specify is the poloidal current Unlike fi,p(\l/), which was determined from the 

frozen-in condition on the disk surface, the function cannot be explicitly prescribed 

as an arbitrary function on any given surface. Instead, it must be somehow determined 
self-consistently together with the solution \l/(r, 9) itself. This means that there must be one 
more condition that we have not yet used. And indeed, this additional condition is readily 
found — it is the (inner) light-cylinder regularity condition. Let us look at it more closely. 

As can be easily seen from the Grad-Shafranov equation (27), the light cylinder, defined as 
the surface where 

D = =>- a = a LC = \6£l\& , (46) 

is a singular surface, because the coefficients in front of both the r- and 0- second-order 
derivatives of \& vanish there. Physically speaking, the light cylinder is the surface where the 
locally-measured rotational velocity of the magnetic field lines with respect to the ZAMOs 
is equal to the speed of light, v-q^ = c, and where E = B poi in the ZAMO frame. In general 
relativity there are two light cylinders, the inner one and the outer one. The outer light 
cylinder is just a direct analog of the pulsar light cylinder; it is crossed by rotating field lines 
that are open and extend to infinity. In our problem, we are interested in the closed field 
lines, i.e., those reaching the event horizon. These field lines cross the so-called inner light 
cylinder, whose existence is a purely general-relativistic effect, first noticed by Znajek (1977) 
and by BZ77. 

Because the inner light cylinder is a singular surface of equation (27), in general this equation 
admits solutions that are not continuous or continuously different iable at the light cylinder. 
Such solutions, while admissible mathematically, are not physically possible. Thus, we sup- 
plement our mathematical problem by an additional physical requirement that the solution 
be continuous and smooth across the light cylinder surface. In particular, this means that 
the 1st and 2nd derivatives of ^ must be finite there. Correspondingly, one can just drop 
all the terms proportional to D when applying equation (27) at the light cylinder and keep 
only the terms involving the derivatives of D. The result can be formulated as an expression 
that determines the function namely: 

-//'(*) = ^ r (d r D)\ LC + \y e (d e D)\ LC + <M^7 2 (W) 2 | LC , (47) 

where \&, r, and are taken at the light cylinder: 



tt = *Lc(0) = tf[r L c(0),0], 



(48) 
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and the function r-^ciO) — the shape of the light cylinder surface — is determined implicitly 
by equation (46). This approach was first used successfully at the outer light cylinder by 
Contopoulos et al. (1999) in the context of pulsar magnetspheres. In the black hole problem, 
it was first used by Uzdensky (2004) for the Schwarzschild case. 

Let us now discuss how one can use the light cylinder regularity condition (47) to determine 
in practice. Conceptually, one can think of this condition as follows. Suppose one 
starts by fixing all the other boundary and regularity conditions in the problem [including 
the choice of f2p(\l/)]. Then, for an arbitrarily chosen function one can regard the 

condition (47) as a mixed-type, Dirichlet-Neumann boundary condition because it can be 
viewed as a quadratic algebraic equation for, say, the first radial derivative. Thus, if 
is given, one can express ^ r \hc i n terms of ^lc an d ^e\hc- Next, one applies this condition 
separately on each side of the light cylinder and gets a complete, well-defined problem in 
each of the two regions separated by the light cylinder. Then, one can obtain a solution in 
each of these regions. Because of the use of the regularity condition (47), each of the two 
solutions is going to be regular near the light cylinder. In general, however, these solutions 
are not going to match each other at r = tlc(^) and the mismatch A^lcW wm depend on 
the original choice of the function This observation suggests a method for selecting a 

unique function one can devise a procedure in which one iterates with respect to 1(9) 

until A^lc becomes zero. The corresponding function is then declared the correct one: 
only with this choice of the solution \l/(r, 9) passes smoothly through the light cylinder. 

The above method for determining is conceptually illuminating and can be easily 

implemented in simple cases. For example, in the case of a uniformly-rotating pulsar magne- 
tosphere, two important simplifications take place. First, the location of the light cylinder is 
known a priori, r^c(0) = c/Q = const, and hence one can choose a computational grid that 
is most suitable for dealing with the light cylinder (e.g., cylindrical polar coordinates with 
some gridpoints lying on the cylinder). Second, because Qp = const, the terms quadratic 
in the derivatives of \l/ disappear, and the task of resolving equation (47) with respect to 
the derivative normal to the light cylinder becomes trivial. These simplifications make the 
procedure described above very practical and it was in fact used successfully by Contopou- 
los et al. (1999) (and repeated later by Ogura & Kojima 2003) to obtain a unique solution 
for an axisymmetric pulsar magnetosphere that was smooth across the outer light cylinder. 

In the problem considered in this paper, however, the situation is much more complicated. In 
particular, the light cylinder's position and shape are not known a priori; instead, they need 
to be determined self-consistently as part of the solution. Also, equation (47) is, in general, 
quadratic with respect to d r ^, and hence one has to deal with the problem of the existence 
of its solutions and with the task of selecting only one of them. Because of this overall 
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complexity, we decided against using this procedure in our calculations. Instead, we chose a 
much simpler and more straight-forward method: we used equation (47) to determine 
[or, rather, the combination Il'fif) that is actually needed for further computations] directly, 
by explicitly interpolating all the terms on the right-hand- side of equation (47). We will 
describe this in more detail in the next section. 



4.3. Numerical procedure 

We performed our calculations in the domain {r e [r#, oo], 9 G [0, vr/2]} on a grid that was 
uniform in 9 and in the variable x = \Jrujr (which enabled us to extend the computational 
domain to infinity). The highest resolution used was 60 gridzones in the ^-direction and 200 
gridzones in the radial (x) direction. To solve the elliptic Grad-Shafranov equation (27), 
we employed a relaxation procedure similar to the one employed by Uzdensky et al. (2002). 
In this procedure, we introduced artificial time variable t and evolved the flux function 
according to the parabolic equation 

™=±f(r,6)(LHS-RHS), (49) 

where LHS and RHS are the left- and right-hand sides of the Grad-Shafranov equation (27), 
respectively, and the factor /(r, 9) was an artificial multiplier introduced in order to accelerate 
convergence in regions where the diffusion coefficients in the x and 9 directions are small 
(e.g., very far away or very close to the horizon). The sign in front of f(r,9) was chosen 
according to the sign of the diffusion coefficient in equation (27): it was plus outside the 
light cylinder (where D > 0) and minus inside (where D < 0). It is clear that any steady- 
state configuration achieved as a result of this evolution is a solution of the Grad-Shafranov 
equation (27). 

Here we would like to draw attention to the following non-trivial problem. During the 
relaxational evolution described by equation (49), the light cylinder generally moves across 
the grid and, from time to time inevitably gets close to some of the gridpoints. This leads to 
the danger, first noted by Macdonald (1984), that some gridpoints will oscillate between the 
two sides of the light cylinder. Indeed, suppose that a given gridpoint P is initially on the 
outer side of the light cylinder (Dp > 0), and so d^p/dt is determined by equation (49) with 
the plus sign. Let us suppose that the resulting evolution of ^ p is such that D P decreases. 
The, after some time one may find that D P has become negative; correspondingly, at the next 
timestep one uses equation (49) with the minus sign and so \&p starts to evolve in the opposite 
direction. Because the value of D at a fixed spatial point P is, locally, a smooth monotonic 
function of ^p, it now starts to increase and may become positive again in one or two 
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timesteps. This leads to rapid small-amplitude oscillations of the light cylinder around some 
gridpoints, instead of a smooth large-scale motion associated with the iteration process. As 
a result, the light cylinder gets "stuck" on these gridpoints and the function rLc(#) becomes 
a series of steps and plateaus instead of a smooth curve. A simple and efficient way to avoid 
this problem turned out to be to update the function D(r,8) not at every timestep but 
rather very infrequently. Although it caused some delay in the convergence of the relaxation 
process, this modification has worked very well in practice, enabling the light cylinder to 
move freely across the grid and to achieve its ultimate smooth shape. 

To implement our relaxation procedure numerically, we used an explicit finite-difference 
scheme with lst-order accurate time derivative and centered 2nd-order accurate spatial 
derivatives. It is also worth mentioning that writing out the full-derivative terms such as 
d r (^ r DA/ p 2 ) as ^ r d r (DA/p 2 ) + (DA/ p 2 )^ rr , etc., and then evaluating them on the grid 
actually worked better than evaluating these full derivatives directly as they are. The initial 
condition — the starting point of our relaxation process — was prescribed explicitly as 

tt(f = 0, r, 9) = ^ s + [tt d (r) - (1 - cos 0) . (50) 

Also, we found it useful to use cubic-spline interpolation of functions I{^>) and f2p(\l/) to 
avoid some small-scale rapid oscillations of the solution. 

Finally, let us describe the particular numerical implementation of the procedure that was 
used to determine the poloidal current J(\P) in our code. As we mentioned at the end of 
the previous section, we used equation (47) explicitly to determine II'(^/) by interpolating 
all the terms on the right-hand side of that equation at the light cylinder. Because the 
light cylinder surface is roughly spherical, it was convenient to represent by a tabular 

function specified on a one-dimensional array {^lc} of the values of \&lc a t the radial rays 
6 = 0i = jhg, where hg is the grid-spacing in 9. Along each of these rays, we first had to 
locate the pair of radially-adjacent gridpoints between which the light cylinder lay. Then we 
used an interpolation of D(r,9) to determine the position r = rLc(^) of the light cylinder 
more precisely and to obtain ty 3 LC = ^lc(0 j '), as we U as the values of the derivatives \I/ r , 
^/g, D r , and D e at the light cylinder for each of the rays. Finally, we used (47) to compute 
the value of II'(^>) at each ^ J LC - This is actually not as trivial as it may seem, because the 
condition (47) in such an approach was enforced at all times during the relaxation procedure 
that determined \l/(r, 0), whereas the Grad-Shafranov equation itself was satisfied only after 
convergence had been reached. Therefore, one had to exercise extra care, for example, in 
deciding how often needs to be updated. We found that it was necessary to update 
J(\l/) only fairly infrequently during our relaxation procedure. 



4.4. Results 



The single most important result of the present study is presented in Figure 4. This figure 
shows where in the two-dimensional (a, \& s ) parameter space force-free solutions exist and 
where they do not. Filled circles on this plot represent the runs in which convergence was 
achieved (allowed region), whereas open circles correspond to the runs that failed to converge 
to a suitable solution (forbidden region). The boundary a max (^ s ) between the allowed and 
forbidden regions is located somewhere inside the narrow hatched band that runs from the 
lower left to the upper right of the Figure [the finite width of the band represents the 
uncertainty in a max (^ r s ) due to a limited number of runs]. As we can see, a max (^ s ) is a 
monotonically increasing function. In particular, in the limit \I/ S — > 0, a max indeed scales 
linearly with and hence is inversely proportional to R s = r m ^/ t ot/^s, in full agreement 
with our expectations presented in § 3. However, this linear dependence no longer holds for 
finite values of ty s (and hence of a max ). 

In order to study the effect that black hole spin has on the solutions, we concentrate on several 
values of a for a fixed value of \I/ S . In particular, we choose \I/ S = 0.5 [that corresponds to 
R s = 2r in (a)] and considered four values of'a: a = 0, a = 0.25, a = 0.5, and a = 0.7. Thus, 
Figure 5 shows the contour plots of the poloidal magnetic flux for these four cases. We 
see that the flux surfaces inflate somewhat with increased a, but this expansion is not very 
dramatic, even in the case a = 0.7, which is very close to the critical value a max (\I / s = 0.5) 
that corresponds to a sudden loss of equilibrium. We note that this finding is completely in 
line with our discussion in § 3. 

The next three Figures present the plots of three important functions that characterize the 
solutions. In each Figure there are four curves corresponding to our selected values a = 0, 
0.25, 0.5, and 0.7 for = 0.5. Figure 6 shows the event-horizon flux distribution \l/o(#); 
Figure 7 shows the poloidal-current function and Figure 8 shows the position of the 

inner light cylinder described in terms of the lapse function «lc(^)- 

In Figure 6 we also plot ^5/^(6) corresponding to the simple split-monopole solution with 
uniform radial field at the horizon. Note that on the horizon we have A = and hence 
£ = r 2 H + a 2 = 2Mr H ; therefore 

Bf = —* = —}—*g = —± ^^e, r = r H . (51) 

wp £sm0 2Mr H smO 

Thus, Bf(9) = const corresponds to ^q\6) = ^ s + (^tot — ^sXl — cos#), independent of a. 
This function is plotted in Figure 6 (the dashed line) for comparison with the actual solutions. 
We see that the deviation from ^^(8) becomes noticeable only when a approaches 1. 
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Figure 8 shows chlc(9). An interesting feature here is that the light cylinder reaches the event 
horizon at some intermediate angle < 9 co < n/2 for small values of a. This is because, 
when a < 0.359.. M, the inner edge of a Keplerian disk rotates faster than the black hole; 
correspondingly, somewhere in the disk there exists a corotation point r co > r in such that 
^A"(^co) = The field line \l/ co threading the disk at this point corotates with the black 
hole. Therefore, at the point 9 co where this line intersects the horizon, we have 8Q — 0, and 
so this point (r = r H , 9 = 9 co ) has to lie on the light cylinder. The location 9 co of this point 
moves towards the equator when a is increased and reaches it at a = 0.359.. M. For larger 
values of a, the entire disk outside of the ISCO rotates slower than the black hole and the 
light cylinder touches the horizon only at the pole = 0. 

Finally, we also computed all the electric and magnetic field components and checked that 
E 2 < B 2 everywhere outside the horizon. 



5. Astrophysical Implications/Consequences 

In this section we'll discuss the exchange of energy and angular momentum between the black 
hole and the disk. Apart from the question of existence of solutions, this issue is one of the 
most important for actual astrophysical applications. Fortunately, once a particular solution 
describing the magnetosphere is obtained, computing the energy and angular momentum 
transported by the magnetic field becomes very simple. 

Indeed, according to MT82, angular momentum and red-shifted energy (i.e., "energy at 
infinity" ) are transported along the poloidal field lines through the force- free magnetosphere 
without losses. Thus, the amount of angular momentum AL transported out in a unit of 
global time t through a region between two neighboring poloidal flux surfaces, \P and ^+A^, 
as given by equation (7.6) of MT82 (modified to suit our choice of notation), is 

and the red-shifted power — flux of redshifted energy per unit global time t — is expressed as 

AP = --n F I A^ , (53) 
2 

(see eq. [7.8] of MT82). 



Then, taking into account the contributions from both hemispheres and both sides of the 
disk, we can compute the total magnetic torque exerted by the hole onto the disk per unit t 
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I(V)dV , (54) 



as 

I'tot 

dL 

~dt ~ ~ 

and, correspondingly, the total red-shifted power transferred from the hole onto the disk via 
Poynting flux is 

*tot 

P = - J ft F (tf)J(tf)dtt. (55) 

Next, since in our problem we have an explicit mapping (42) between \1> and the radial coor- 
dinate r on the disk, we can immediately write down expressions for the radial distributions 
of angular momentum and red-shifted energy deposited on the disk per unit global time: 

dAL(r) dV d 

— — — = -I[V d (r)\ — dr , (56) 

and 

AP(r) = -Q K (r) J[tf d (r)] ^ dr . (57) 

Figures 9 and 10 show these distributions for our selected cases a = 0.25, 0.5, and 0.7 for 
fixed ty s = 0.5. We see that in the case a = 0.25 there is a corotation point r co on the disk 
such that fidisk > inside r co and fidisk < &h outside r co . Correspondingly, both angular 
memontum and red-shifted energy flow from the inner (r < r co ) part of the disk to the black 
hole and from the hole to the outer (r > r co ) part of the disk. At larger values of a, however, 
the Keplerian angular velocity at r = r in is smaller than the black hole's rotation rate and 
there is no corotation point; correspondingly, both angular momentum and redshifted energy 
flow from the hole to the disk. Also, as can be seen in Figures 9 and 10, the deposition of 
these quantities becomes strongly concentrated near the disk's edge, especially at higher 
values of a. 

Next, Figures 11 and 12 demonstrate the dependence of the total integrated angular mo- 
mentum and red-shifted energy fluxes (54)-(55) on the black hole spin a for a fixed value of 
= 0.5. We see that both quantities are negative at small values of a (meaning a transfer 
from the disk to the hole,) but then increase and become positive at larger a. The angular 
momentum transfer rate depends roughly linearly on a, whereas the red-shifted power P(a) 
grows even faster, especially at large values of a. It is also interesting to note that the two 
quantities go through zero at slightly different values of the spin parameter: dL/dt becomes 
zero at a 0.23, while P = at a as 0.26. This means that it possible to have the total 
angular momentum flow from the hole to the disk and the total power flow from the disk to 
the hole at the same time. 
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6. Conclusions 

In this paper we investigated the structure and the conditions for the existence of a force-free 
magnetosphere linking a rotating Kerr black hole to its accretion disk. We assumed that 
the magnetosphere is stationary, axisymmetric, and degenerate and that the disk is thin, 
ideally conducting, and Keplerian and that it is truncated at the Innermost Stable Circular 
Orbit. Our main goal was to determine under which conditions a force-free magnetic field 
can connect the hole directly to the disk and how the black hole rotation limits the radial 
extent of such a link on the disk surface. 

We first introduced (in § 3) a very simple but robust physical argument that shows that, 
generally, magnetic field lines connecting the polar region of a spinning black hole to arbi- 
trarily remote regions of the disk cannot be in a force-free equilibrium. The basic reason for 
this can be described as follows. Since the field lines threading the horizon have to first cross 
the inner light cylinder, and since they generally rotate at a rate that is different from the 
rotation rate of the black hole, these field lines have to be bent somewhat. In other words, 
they develop a toroidal magnetic field component, just like the open field lines crossing the 
outer light cylinder in a pulsar magnetosphere. In the language of the Membrane Paradigm 
(see Thorne et al. 1986), this toroidal field is needed so that the field lines could slip re- 
sistively across the stretched event horizon. The next step in our argument is to look at 
those field lines that connect the polar region of the horizon to the disk somewhere far away 
from the black hole. In a force-free magnetosphere, toroidal flux spreads along field lines to 
keep the poloidal current / ~ B^azu constant along the field. Then one can show that the 
outward pressure of the toroidal field generated due to the black hole rotation turns out to 
be so large that it cannot be confined by the poloidal field tension at large enough distances. 
In other words, the field lines under consideration cannot be in a force-free equilibrium. 
Furthermore, one can generalize this argument to the case of closed magnetospheres of finite 
size and derive a conjecture that the maximal radial extent R max of the magnetically-coupled 
region on the disk surface should scale inversely with the black hole spin parameter a in the 
limit a — > 0. 

In order to verify this hypothesis and to study the detailed structure of magnetically-coupled 
disk-hole configurations, we have obtained numerical solutions of the general-relativistic 
force-free Grad-Shafranov equation corresponding to partially-closed field configurations 
(shown in Fig. 3). This is an nonlinear 2nd-order partial differential equation for the poloidal 
flux function \l/(r, 6) and it is the main equation governing the system's behavior. 

An additional complication in this problem arises from the need to specify two free functions 
that enter the force-free Grad-Shafranov equation; these are the field-line angular velocity 
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^f(^) and the poloidal current Because all the field lines are assumed to be frozen into 

the disk, the first of this functions is determined in a fairly straightforward way. Namely, for 
any given field line Qp^) is just the Keplerian angular velocity at this line's footpoint 
on the disk. Specifying the poloidal current, on the other hand, is a much more difficult 
and nontrivial task. The reason for this is that it cannot be just prescribed explicitly on 
any given surface and one should look more thoroughly into the mathematical nature of the 
Grad-Shafranov equation itself to determine In particular, the most important feature 

of the Grad-Shafranov equation in this regard is that it becomes singular on two surfaces, 
the event horizon and the inner light cylinder. This observation is very useful because one 
can impose a physically-motivated regularity condition at each of these surfaces. One of 
the most important ideas in our analysis is that one can use the light- cylinder regularity 
condition to determine, using an iterative procedure, the poloidal current similar to 

the way it was done by CKF99 in the context of pulsar magnetospheres. 

As for the singularity at the event horizon, it is also very important. Basically, it tells us that 
it is not possible to prescribe an arbitrary boundary condition at the horizon; instead, one 
can only impose a certain physical condition of regularity there. When combined with the 
Grad-Shafranov equation itself, this regularity requirement results in a single relationship 
(historically known as the horizon boundary condition, first derived by Znajek 1977) between 
three functions: the horizon flux distribution ty (9), an d the two free functions, Qf(^) 
and 1(9) (e.g., Beskin 1997; Beskin & Kuznetsova 2000). What's important is that there 
are no other independent relationships that can be specified on this surface. In practical 
terms, this means that this condition should be used to determine the function 9 (9) in 
terms of ^f(^) and 1(9), which therefore must be determined outside the horizon. This 
fact helps to alleviate some of the causality issues raised by Punsly (1989, 2001, 2003) and 
by Punsly & Coroniti (1990). 

Since one of the goals of this work was to study the dependence -R max (a), we performed a 
series of computations corresponding to various values of two parameters: the black hole 
spin parameter a and the the magnetic link's radial extent R s on the disk surface (the 
field lines anchored to the disk beyond R s were taken to be open and non-rotating). At 
the same time, the disk boundary conditions were kept the same in all these runs, namely, 
^d{ r ) — ^ / tot r m( a )/ r - Therefore, varying the value of R s for fixed a was equivalent to varying 
the amount 9 S of open magnetic flux threading the disk. 

Whereas for some pairs of values of a and 9 S we were able to achieve a convergent force-free 
solution, for others we were not. Thus, as one of the main results of our computations, 
we were able to chart out the allowed and the forbidden domains in the two-parameter 
space (a, ^ s ). The boundary between these two domains is a curve a mSuX (9 s ), which can be 
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easily remapped into the curve R max (a). As can be seen in Figure 4, this is a monotonically 
rising curve with the asymptotic behavior a max oc *ff s as — > 0, which is in line with our 
predictions. 

We also computed the total angular momentum and red-shifted energy exchanged in a unit 
of global time t between the hole and the disk through magnetic coupling. We studied 
the dependence of these quantities on the black hole spin parameter a and found that the 
angular momentum transfer rate rises roughly linearly with a; it is negative for small a 
(meaning the angular momentum transfer to the hole) and reverses sign around a « 0.23 
(for \l/ s = 0.5^tot)- The total energy transfer increases with a at an accelerated (i.e., faster 
than linear) rate, especially at larger values of a; it is also negative at small a, but becomes 
positive around a = 0.26. This means that there is a narrow range 0.23.. < a < 0.26.. where 
the integrated angular momentum flows from the hole to the disk, whereas the integrated 
red-shifted energy flows in the opposite direction. 

Finally, we note that, in the case of open or partially-open field configuration responsible for 
the Blandford-Znajek process, one has to consider magnetic field lines that extend from the 
event horizon out to infinity. Since these field lines are not attached to a heavy infinitely 
conducting disk, their angular velocity J1f(\I>) cannot be explicitly prescribed; it becomes 
just as undetermined as the poloidal current they carry. Fortunately, however, these 
field lines now have to cross two light cylinders (the inner one and the outer one). Since each 
of these is a singular surface of the Grad-Shafranov equation, one can impose corresponding 
regularity conditions on these two surfaces. Thus, we propose that one should be able 
to devise an iterative scheme that uses the two light-cylinder regularity conditions in a 
coordinated manner to determine the two free functions fi(^) and I(^f) simultaneously, as a 
part of the overall solution process. At the same time, the regularity conditions at the event 
horizon and at infinity could be used to obtain the asymptotic poloidal flux distributions 
at r = r H and at r — > oo, respectively. We realise of course that iterating with respect to 
two functions simultaneously may be a very difficult task. This purely technical obstacle (in 
addition to having to deal with the separatrix between the open- and closed-field regions) is 
the primary reason why, in this paper, we have restricted ourselves to a configuration which 
has no open field lines extending from the black hole to infinity. We leave this problem as a 
topic for future research. 

It is possible that, instead of solving the Grad-Shafranov equation itself, the easiest and 
most practical way to achieve a stationary solution will be to use a time- dependent rela- 
tivistic force-free code, such as one of those being developed now (Komissarov 2001, 2002a, 
2004a; MacFadyen & Blandford 2003; Spitkovsky 2004; Krasnopolsky 2004, private commu- 
nication). 
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i \ disk 
black hole 

Fig. 1. — Schematic drawing of an open black hole - disk magnetosphere, commonly associ- 
ated with the Blandford-Znajek (BZ77) process. Rotational energy and angular momentum 
are extracted from both the black hole and the disk and are transported away by the magnetic 
field. 
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black hole 

Fig. 2. — Schematic drawing of a fully-closed black hole - disk magnetosphere. Energy and 
angular momentum are exchanged between the hole and the disk through a direct magnetic 
link. 
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Fig. 3. — Schematic drawing of a black hole - disk magnetosphere with a radially-limited 
magnetic connection. Here, only the inner part of the disk is coupled magnetically to the 
hole (closed field region), whereas the field lines attached to the outer part of the disk are 
open and extend to infinity. 
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Fig. 4. — The 2-D parameter space (a, \& s ). Filled (black) circles correspond to the runs in 
which a stationary force-free solution has been obtained, while open (white) circles corre- 
spond to the runs that failed to converge to a stationary solution. The shaded band running 
diagonally across the plot represents the function a max (i? s ), the maximal value of a for which 
a force-free magnetic link can extend up to a given radial distance R s on the disk. 
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open total 

Fig. 5. — Contour plots of the magnetic flux function \I/(r, 9) for four values of the black hole 
specific angular momentum: a = 0.0, 0.25, 0.5, and 0.7; the amount of open poloidal flux in 
all cases is \l/ s = 0.5^ t ot (corresponding to R s = 2r in ). 
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Fig. 6. — The horizon poloidal flux distribution ^o(0) for several values of a and a single 
value \l/ s = 0.5\l/ tot of the open magnetic flux (corresponding to R s = 2r in (a). 
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Fig. 7. — Poloidal current / as a function of poloidal magnetic flux ^ for several values of a 
and a single value \l/ s = 0.5\l/ tot . 
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Fig. 8. — The position of the inner light cylinder represented by the lapse function a\ J c{9) 
for several values of a. The light cylinder touches the horizon at the pole 6 = and, for 
a < 0.36, at the point 9 = 9 co where the corotation field line ^/ co = ^d(r co ) intersects the 
horizon. 




Fig. 9. — Radial distribution d 2 L/drdt of the magnetic torque per unit radius r on the disk 
surface for a = 0.25, 0.5, and 0.7. 




Fig. 10. — Radial distribution dP/dr of the red-shifted power per unit radius r on the disk 
surface for a = 0.25, 0.5, and 0.7. 
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Fig. 11. — The dependence of the total magnetic torque between the black hole and the disk 
on the hole's spin parameter a for fixed R s = 2r in . 
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Fig. 12. — The dependence of the total red-shifted power exchanged magnetically between 
the black hole and the disk on the hole's spin parameter a for fixed R s = 2r in . 



